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ABSTRACT 

We develop a ID solar- wind model that includes separate energy equations for the electrons and protons, pro- 
ton temperature anisotropy, coUisional and collisionless heat flux, and an analytical treatment of low-frequency, 
reflection-driven, Alfven-wave (AW) turbulence. To partition the turbulent heating between electron heating, 
parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave 
damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increas- 
ing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for 
either instability. We numerically integrate the equations of the model forward in time until a steady state 
is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of 
observations, supporting the idea that AW turbulence plays an important role in the origin of the solar wind. 
Subject headings: solar wind — Sun: corona — turbulence — waves 
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1. INTRODUCTION 

The first theory for the origin of the solar wind was de- 
yeloped by E. N. Parker in the 1950s and 1960s. Parked 
(119581 11965 ^ based his analysis on a single-fluid, hydrody- 
namic model with electron thermal conduction. Although he 
obtained smooth, transonic solutions in agreement with some 
solar-wind observations, his work was unable to account for 
the large flow velocities and proton temperatures measured in 
fast-solar-wind streams near Earth. Subsequent studies gen- 
eralized the Parker model to include separate energy equa- 
tions for the protons and electrons ( Hartle & Sturrock 1968), 
temperature anisotropy ( Leer &_Axford 1972; Whang 1972), 
super-radial expansion of the magnetic field (Holzer & Leer 
[1980), colUsionless heat flux (lHollweail974i il976i) . and en- 
ergy and mome ntum deposition by Alfven waves (AWs) and 
AW t urbulence ( Alazraki & Couturier 1 971l:lBelcher & DavisI 
[MIlLHollweg 1973a b; Tu 1987, 1988). 

The idea that AW turbulence plays an important role in the 
origin of the solar wind was originally proposed bv lColernaiil 
(0968) and is consistent with a number of observations. For 
example, data obtained from the Solar Optical Telescope on 
the Hinode spacecraft reveal the presence of ubiquitous, AW- 
like motions in the low corona carrying an energy fl ux suffi- 
cient to power the solar wind ( De Pontieu et alj|2007 h. Perva- 
sive AW-like fluctuations are also seen at higher altitudes in 
the corona in observations from the Coronal Multichan nel Po- 
larimeter of the National Solar Observatory (Tomczv k et alj 
|2007). Voyager, Helios, Wind, and other spacecraft have 
measured broad-spectrum fluctuations in the magnetic field, 
electric field, and flow velocity in the solar wind, demon- 
strating that AW t urbulence is present throughout the inter- 
planetary medium jTu & Marschlll995l:lGoldsteiii et al.ll 19951: 
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iBale et aO2005h . Moreover, the amplitudes of these turbulent 
fluctuations are sufficient to explain the heating rates that have 
been inferred from the proton and electron temperatur e pro- 
files ( ISmith et al. 2001; MacBride et al. 2008; Cranme retaLl 
2009; Staw arz et al.l i2009). Radio-scintillation observations 
of solar-wind density fluctuations place upper limits on the 
turbulent heating rate in the solar wind that are consistent with 
solar- wind heating by AW turbulence (Harmon & Coles 2005; 
Chandran et al 2009; but see also Spangler 2002). In addi- 
tion, Faraday rotation of radio transmissions from the Helios 
spacecraft show that magnetic-field fluctuations in the corona 
at heliocentric distances between 2Rq and 15Rq are consis- 
tent with models in which the solar wind is driven by AW 
turbulence dHoUweg et al.|[T982ll20Toh . 

However, it is not clear that heating by AW turbulence can 
explain observations of ion temperature anisotropics. Mea- 
surements from the Helios spacecraft have shown that T±p > 
7j|p in the core of the proton velocity distribution in low-p||p 
fast-wind streams, where 



87tnA:B7||] 

B2 



(1) 



n is the proton density, feB7||p/2 is the average energy per pro- 
ton in thermal motions parallel to the magnetic field B, and 
^b7±p is the average ene rgy per p roton in thermal motions 
perpendicular to B ( Marsch et alj [2004). Similarly, remote 
observations from the Ultraviolet Coronagraph Spectrome- 
ter (UVCS) on the Solar and Heliospheric Observatory have 
shown that Tj_ ^ 7|| for O"*"^ ions in coronal holes dKohl et alj 
Il998t iLi et al 1998': 'Ant onucci et a l."2000^. These observa- 
tions pose a challenge to solar-wind models based on AW tur- 
bulence because the AW energy cascade is anisotropic, trans- 
porting AW energy primarily to small scales measured per- 
pendicular to B and only weakly to small scales measured 
parallel to B (Sheba lin et alJll983h . As a consequence, very 
little energy cascades to high frequencies comparable to the 
ion cyclotron frequencies at wh ich waves can dissipate via 
resonant cyclotron interactions dOuataertl 119981) . which are 
the only route to perpendicular ion heating in a colli sion- 
less p lasma within the framework of quasilinear theory dStixl 
fT992h . 
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A number of studies have gone beyond quasilin- 
ear theory to show that low-frequency AW turbulence 
can lead to perpendicular ion heating e ven in the ab- 
sence of a^ cyclo t ron resonance (McChes nev et alJ 119871: 



Chenetal. '200 r, 'Johnson & Cheng 2001; DmitruketaJJ 
2004: Voi tenko & Goossens 2004: Bourouaine et al. 2001 
Parashar et al. 2009: Chandran et al. 2010: Chandran 201(|l. 
In this paper, we incorpo rate results from one of these 
studies jChandran et alj |2010i) into a quantitative solar-wind 
model to investigate the extent to which low-frequency AW 
turbulence can explain the observations of anisotropic pro- 
ton temperatures discussed above. The model we have devel- 
oped also includes the non-WKB reflection of AWs, proton 
and electron heat flux in both the collisional and collisionless 
regimes, and enhanced pitch-angle scattering when the proton 
temperature anisotropy is sufficiently large to excite mirror, 
cyclotron, or firehose instabilities. We describe this model in 
detail in Section |2] and discuss the numerical method we use 
to solve the equations of the model in Section|3] In Section|4] 
we present and analyze a steady-state, fast-wind-like solution 
and compare this solution to a number of observations. In 
Section |5] we present a second steady-state solution and dis- 
cuss our results. We conclude in Section|6]and summarize the 
derivation of the equations of our model in the appendix. 

2. MATHEMATICAL FRAMEWORK 

In the solar wind, the Debye length and proton gyrora- 
dius Pp are vastly smaller than the length scales over which 
the bulk solar-wind properties vary appreciably. Because of 
this, the large-scale structure of solar-wind plasma can be 
rigorously described by taking the limit of the Vlasov and 
Maxwell equations in which Pp — > and — > 0. This limit, 
described in more detail in the appendix, is some times re- 
ferred to as K ulsrud's collisionless magnetohydrodynamics 
(lKulsrudlll983"h . This name is somewhat misleading, in that 
the resulting description is still kinetic in nature. In particu- 
lar, two of the variables in Kulsrud's theory are the reduced 
proton and electron distribution functions, /p and /e, which 
are independent of the gyrophase angle in velocity space and 
evolve in time according to the guiding-center Vlasov equa- 
tion (Equation ( IA9I )). 

We take Kulsrud's collisionless magnetohydrodynamics 
(MHD) as the starting point for our study, including the sim- 
plifying assumption that the p lasma consis t s of ju st protons 
and electrons. We then follow ISnyder et aH (11997 ) in adding 
a collision operator to the guiding-center Vlasov equation and 
taking velocity moments of this equation to obtain a hierarchy 
of fluid equations, as described in more detail in the appendix. 
To close this set of fluid equations for the protons, we set /p 
equal to a bi-Maxwellian when evaluating the various fourth 
moments of /p. This approach offers a simple and, we expect, 
reasonably accurate way to solve for the proton heat flux in 
the presence of temperature anisotropy in both collisional and 
collisionless conditions. The same general equations (with 
differing treatments of the collision ter ms) were derived in 
different ways by Endeve & Leer (2001'). lLie-Svendsen et alj 
((2001), and Ramos (20031) . For the electrons, we adopt the 
simplifying assumption that /e is isotropic in velocity space 
and close the fluid equations by specifying the electron heat 
flux in terms of lower-order moments of /e. The most general 
versions of the resulting equations are given in the appendix. 
In this section, we specialize these equations to our ID solar- 
wind model and add terms to incorporate heating by AW wave 



turbulence, acceleration by the AW pressure force, and tem- 
perature isotropization by firehose and mirror instabilities. 

To simplify the analysis, we neglect the Sun's rotation, take 
the background magnetic field to be fixed, and solve the fluid 
equations within a narrow, open magnetic flux tube centered 
on a r adial magnetic field line. We follow Kopp & Holze^ 
(119761) in taking the cross sectional area of this flux tube to be 



where r is heliocentric distance. 
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/i = 1- (/max - l)exp[(/;0-7;i)/o], anda0,/max,^i, ando 
are constants. The function / increases from 1 at r = Rrr, to 
/max at r ^ with most of the variation in / occurring be- 
tween r = Ri ~ a and r = Ri + a. To reduce the number of 
free parameters, we set 



a^Ri 



(4) 



The inner radius of our model corresponds to the coronal base 
just above the transition region. We set this radius equal to Rq, 
neglecting the thickness of the chromosphere and transition 
region. By flux conservation, the strength of the magnetic 
field B satisfies 

B=^, (5) 
a 

where Bq is the magnetic field strength at the coronal base. 
We take the cross sectional area at the coronal base, Op,, to 
be <C R% and <C R^, so that the flux tube is thin. As a con- 
sequence, the magnetic field direction is approximately radial 
everywhere within the flux tube. We thus set 
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where we have taken the magnetic field to be pointing away 
from the Sun. The condition V • -B = then gives 



_ - Ida 
Vb^- — . 
a or 



(7) 



We take the solar- wind outflow velocity to be everywhere par- 
allel to the magnetic field, 

U = Ub, (8) 
and define the Lagrangian time derivative 

d d d 
dt-di^^Tr- 

There are eight dependent variables in our model: the pro- 
ton (or electron) number density n, the proton (or electron) 
outflow velocity U, the electron temperature Tf., the perpen- 
dicular and parallel proton temperatures T^p and 7||p, the pro- 
ton heat fluxes q±p and <7||p, and the energy density of AWs 
propagating away from the Sun £w The proton heat flux ^xp 
is a flow along B of perpendicular proton kinetic energy: no 
heat flows across the magnetic field in the model. The eight 
dependent variables of the model depend upon time t and 
a single spatial coordinate r and satisfy the following eight 
equations; 

dn n d , 

-r^--^{aU), (10) 
dt a or 
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Vpn^B(71p - 7||p) +VpenfcB(7; - 7||p), (14) 



^2d_ /qj_p\ _ _ w^b^Ip 371p w^^XpC^Xp ^ ^|p) 3a 
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and dPewdl 19701) 



3£w 
where 



-|-[a(f/ + VA)£w] + f^|-(af/) = -e, (17) 
a or la dr 
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is the Alfven speed, p is the mass density, and Mq is the mass 
of the Sun. Since we are treating the solar wind as a proton- 
electron plasma, we ignore the contribution of alpha particles 
and other particle species to p and set 



P = OTpH. 



(19) 



The quantities Q^, Qxp. and 2||p are, respectively, the elec- 
tron heating rate, the perpendicular proton heating rate, and 
the parallel proton heating rate per unit volume from the dis- 
sipation of AW turbulence (Section lZ4l i. 

G = Ge + G±p + e||p (20) 

is the total turbulent heating rate per unit volume, and is 
the electron heat flux (Section lOT l. The quantity 
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(21) 



is the Coulomb collision frequency for energy exchange be- 
tween protons and electrons (,Schunk.l975i) . where and me 



are the proton and electron masses and InA is the Coulomb 
logarithm, which we take to be 23. We set 

(22) 
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where 



^pp.c 
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(23) 



is the proton-proton Coulomb collision frequency (ISchunkI 
Vinst is a scattering rate associated with small-scale 
plasma waves that are excited when the prot on te mperature 
anisotropy becomes sufficiently large (Section |272] i, and 
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J-p- 



-. (24) 

Equation ( fTOl l expresses the conservation of mass in our ID 
model. Equation (fTTT i is the same as the momentum equa- 
tion in Kulsrud's collisionless MHD (Equation (IA2I) ). ex- 
cept that we have add ed the gravitational acceleration and the 
AW pressure force ( Dewaij[l97 0). In the absence of turbu- 
lent heating, heat flow, and collisions, the right-hand sides 
of Equations ( fT2] i through ( fT4b vanish. In this case, the 
Lagrangian time derivative of the electron specific entropy 
(oc ln(re/n^''^)) vanishes, and the protons obey the double- 
adiabatic theory of Chew, Goldberger, & Low (1956). When 
Vp becomes sufficiently large, 7xp — 7]|p and the collisional 
terms in Equations ( fTSl l and ( fT6b (which are 0= Vp) are much 
larger than the left-hand sides of these equations. In this limit, 
^Xp and <7||p are determined from Equations ( fTSI ) and ( fT6l l 
by balancing the collisional terms against the source terms 
(which contain Txp and/or 7||p but not <7xp or <7||p), and the to- 
tal proton heat flux g'xp + ^||p becomes approximately equal to 
the proton heat flux in collisional transport theory dBraginskiil 
1965). 

Upon multiplying Equation (fTTT i by pt/ and adding the re- 
sulting equation to the sum of Equations (fTSI) . ( fT3t , (fT4l) . 
and (flTb . we obtain a total energy equation. 
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where 



-nkn 



37; 



is the total energy density, and 
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(27) 



is the total energy flux. In steady state, aFtot is independent 
of r, and the total flow of energy into the base of the flux 
tube equals the total flow of energy through the flux-tube cross 
section at all other radii. 

2.1. Electron Heat Flux 

Close to the Sun, n is sufficiently large that the electron 
Coulomb mean free path. 



Wfp 



\Jk^%ln 



(28) 
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is much shorter than the radial distance h — T^/ (dTe/dr) over 
which Tg varies appreciably. The quantity 



Ve =2.9x 10 
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1 cm 



1 eV 



-3/2 
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(29) 



is the electron collision frequency ('Book"1983'). We expect 
the electron heat flux in this nea r-Sun region to be app roxi- 
mately equal to the Spitzer value jSpitzer & Harm|[l953h . 



ge,s 
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where 



KeO = 



1.84 X 10- 
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■ erg s" 



cm 



(30) 



(31) 



Farther from the Sun, A-mfp > It, and Qf. deviates from the 
Spitzer value. We follow Hollweg (1974, 1976) in taking the 
colUsionless heat flux in this region to be approximately 



(32) 



whe re an is a constant that we treat as a free parame- 
ter. iHollwegl {1976) argued that the transition between the 
coUisional and collisionless regimes occurs at the radius at 
which A-mfp — 0.5r. To interpolate smoothly between the two 
regimes, we set the electron heat flux equal to 



where 



(33) 
(34) 



and rH is a constant that we choose to coincide with the radius 
at which A-mfp — 0.5r . For the numerical solutions presented 
in Sectionsl?[and l5.7l we set rw = 5Rq, and confirm post facto 
that A ~ 0.5r at r = rn (see Figure|6]l. 

2.2. Proton Pitch-Angle Scattering from Firehose and Mirror 
Instabilities 



If the proton temperature-anisotropy ratio 



R 



(35) 



becomes either too large or too small, the plasma becomes 
unstable. Spacecraft measurements show that the values of R 
found in the solar wind are bounded from below by the insta- 
bility threshold of the oblique firehose mode an d from above 
by the insta bility thr eshold of the mirror mode (iKasper et al.l 
|2002; Hel Unger et al. 2006; Bale et al. 2009). In particular, 
most of the measured values of R correspond to plasma pa- 
rameters for which Ymax < lO^^Hp, where ymax is the maxi- 
mum growth rate of the oblique firehose or mirror instability 
and iip is the proton cyclotron frequency. The value of R for 
which Ymax — lO^^ilp is approximately 



-0.76 



(36) 



/?m = l+0.77(p||p + 0.016) 
for the mirror instability, and approximately 

/?f-l-1.4(p||p + 0.11)-i (37) 
for the oblique firehose instability dHelUnger et al.ll2006h . 



Presumably, when the plasma becomes unstable, small- 
scale electromagnetic fluctuations grow and enhance the pro- 
ton pitch-angle scattering rate, preventing the temperature 
anisotropy from increasing further. We incorporate this ef- 
fect into our model through the term Vinst in Equation 
with 



Vinst =Voexp 



\2{R-R^) 



Ru 



Voexp 



l2{Rf-R) 



, (38) 



Vo 0.02^JgMq/rI, and Rf = max{R{, IQ-^). A similar 

approach was employed bv lSharma et alj (120061) in numerical 
simulations of accretion flows around black holes. 

2.3. Alfven Wave Turbulence 

The Sun launches different types of waves that propa- 
gate outward into the solar atmosphere. In our model, 
we retain only the non-compressive Alfven wave (AW), 
in part for simplicity and in part because the AW is the 
most promising wave type for transporting en ergy ove r large 
distances into the corona and solar wind ("Barnes' '19661 
Vellietal. 1989; Matthaeus et al. 1999; Suzuki & Inutsuki 
120051: ICranmer & van Ballegooiienil2005 ). For AW fluctua- 
tions, the fluctuating velocity vector 5i> and magnetic field 
vector 5-B lie in the plane perpendicular to Bq. We define the 
Elsasser variables 



(39) 



and, as mentioned previously, take Bq to point away from 
the Sun. In the small-amplitude limit, 2+ fluctuations are 
AWs that propagate with an outward radial velocity ofU + v'a, 
while the fluctuations are AWs that propagate with a radial 
velocity U — va- Near the Sun, U < v'a and fluctuations 
propagate towards smaller r. 

To a good approximation, AWs in the solar corona and 
solar wind can be described within the fr a mework of re- 
duced MHD (Kadomtsev & Pogutse 1974; Str aussI [T976t 
Zank & Matthaeus 1992; Schekochihin et al. 2009^ In re- 
duced MHD, the outward-propagating z+ waves generated 
by the Sun do not interact with one another. However, 
we assume that most of the AW energy is at periods of 
tens of minutes to hours, which makes the wavelengths 
in the radial direction sufficiently long that the AWs un- 
dergo significant non-WKB reflection, converting some of the 
z+ w aves into waves (Heine mann & Olbert 1980c iVelUl 
119931) . Interactions between z+ and fluctuations then 
cause wave energy to cascade fror n large scales to small 
scales (Iroshnikov 1963; Kraichn anlll965t IVelli et all 1 19891; 
iMatthaeus et al. 1999i Ivan Ballegooiien et alJl201 IV At suf- 
ficiently small scales, the z* fluctuations dissipate. Although 
some reflection occurs, we assume that 



z <z^ 



(40) 



and neglect the contribution of z to the wave energy den- 
sity which is then given by 



(41) 



where is the rms amplitude of z^ fluctuations. 

To describe the cascade of wave energy in the pres- 
ence of wave reflections, we adopt the phenomenological 
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model of | Dm itruk e t alj (120021) . which was later extended by 
IChandran& HollweilSoolF to account for the solar-wind 
outflow velocity. The essence of these models is to bal- 
ance the rate at which waves are produced by wave re- 
flections against the rate at which the waves cascade 
and dissipate via interactions with z+ waves. This balance 
lea ds to the following estima te for the rms amplitude of 
(IChandran & Hollwegll2009l) : 



8vA 



(42) 



where Lx is the correlation length (outer scale) of the Alfvenic 
fluctuations in the plane perpendicular to Bq. The rate at 
which energy cascades and dissipates per unit volume is then 



Q = 



CdP^rms (^rms) 

4L, 



(43) 



where is a dimensionless number Since our estimate 
of is proportional to Li_, the value of Q in equation (03]) is 
independent of Lx- Because of Equation ( l40l i. we have omit- 
ted a term ztm& iz^ ^)"^ that is some ti mes included in the 
turbulent heating rate (iHossain et al.lll995.) . 

2.4. Proton and Electron Heating Rates 

In Equation ( [TtI i. the rate Q at which energy is drained 
from the AWs equals the energy cascade rate given in Equa- 
tion (l43l l. which is determined by the "large-scale quantities" 
Grins' ^rms' ^i^^ ^x- ^he AW energy that cascades to 

small scales dissipates, contributing to turbulent heating, but 
the way that Q is apportioned between 2e, 6±p> 2||p de- 
pends upon the mechanisms that dissipate the fluctuations at 
length scales <C ix- In this section, we describe how we di- 
vide the turbulent heating power between Q^, 2xp^ and 2||p 
using results from the theories of linear wave damping and 
nonlinear stochastic heating. 

Nonlinear interactions between counter-propagating AWs 
cause AW energy to cascade primarily to larger A;x and 
only weakly to larger where k^^ and are wavevector 
components perpendicular and parallel to (Shebalin et ajj 
1983; Goldreich & Sridhar 1995; N2 & Bhattachariee 199^ 



Galtier et al.ll2000l) . This cascade does not transfer AW en- 



ergy efficiently to higher frequencies (the AW frequency 
being k^\v/C), and thus cyclotron damping is not an im- 
portant dissipation mechanism for the anisotropi c AW 
cascade (Cranmer & van Ballegooiien 2003; Howe Fet al.l 
|2008a). There may be other mechanisms in the solar wind 
that generate AWs with sufficiently h igh frequencies that 
the AWs undergo cy clotron damping dLeamon et al.l 119981: 
[Hamilton e t al.l l2008h . such as a turbulent cascade involv- 
ing compressive waves dChandran 2005, 2008b; Yoon & FanJ 
|2W)8|) or i nstabilities driven by proton or alpha-particle 
beams (iGomberoff et al...l996: .Helhnger & Travmcek.201 iT) . 
However, we do not account for these possibiUties in our 
model. 

When AW energy cascades to A:xPp — 1, the ca scade tran- 
sition s to a kinetic A lfven wave (KAW) cascade (Bale et al.' 
2005t lHowes etal.1 2008a b; Schekochihin et al. 2009; 
Sahraoui et al ■I I2009h . and the KAW fluct uations undergo 



Landau damping and transit-time dam ping dOuataertl 119981: 
iGruzinovl 119981: iLeamon et al.l 119991) and dissipation via 
stochastic heating (McCh esnev et al.ll 19871 IChen et al.ll200U 
Hohnson & Cheng. .200 1.) . Some of the turbulent energy 



dissipates at A;xPp — 1, and some of the turbulent energy 
cascades to, and then dissipates at, smaller scales. Before 
describing the details of how we incorporate dissipation 
into our model, we first summarize our general approach. 
We make the approximation that the dissipation occurs in 
two distinct wavenumber ranges: A:xPp ^ 1 and A;xPp ^ 1. 
We divide the total dissipation power between these two 
wavenumber ranges by comparing the energy cascade time 
scale and damping time scale at fcxPp = 1 (see Equation ( l53T l 
below). We divide the power that is dissipated at ^xPp ^ 1 
between Q^, Qxp. and 2||p by comparing the damping rates at 
A;xPp = 1 associated with three different dissipation mecha- 
nisms, each of which contributes primarily to either Q^, Qi_^, 
or 2||p- We then assume that all of the power that dissipates 
at A:x ^ Pp ' does so via interactions with electrons, thereby 
contributing to 2e- 

We define Ye and Yp to be the electron and proton contri- 
butions to the linear damping rate of KAWs at A:xPp — 1, 
where Pp is the proton gyroradius. Using a numerical code 
that s olves the full hot-plasma dispersion relation (Ouataerl 
we have calculated Ye and Yp for a range of plasma 
parameters, assuming isotropic proton and electron tempera- 
tures. For 10"^ < Pp < 10, 1 < Tp/T; < 5, and |A:||Va| < iip, 
our results are well approximated by the following formulas: 



I^iivaI 



0.01 



7; 

TpPp 



1/2 



1 



-0.17pi-3 



1 + (2800pe) 



-1.25 



and 



^11 V/ 



= 0.08 



1/4 



po-^exp 



1.3 



(44) 



(45) 



where pp = ?,%nk^T-p / B}^ and Pe = 'Siiink^T^ / Bq. In Figure[Tl 
we compare Equations ( |44| ) and ( 1451) with our numerical solu- 
tions for the case in which Tp — 2Te. 
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Fig. 1. — Solid line gives the e lectr on contribution to the KAW damping 
rate at k±Pp = I from Equation I44t . and the dotted line gives t he p roton 
contribution to the KAW damping rate at Axpp = 1 from Equation i45i . The 
circles and diamonds are, respectively, the electron and proton contributions 
to the KAW damping rate at kj^Pp = 1 in our numerical solutions to the full 
hot-plasma dispersion relation for Maxwellian electrons and protons. 
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At kj^Pp 1, AW/KAW turbulence has a range of 
^11 values. However, we approximate the linear proton 
and electron damping rates by assigning a single effective 
I to the spectrum at k±Pp = 1 g iven by the critical- 
balance conditio n (Higdon 1984; Goldr eich & SridhaiilT99l 
ICho & LazarTanli2004;,Boldyrev.2006) 



l^lllvA^ 



where 



tc = 



P5v| 
Q 



(46) 



(47) 



is the energy cascade time at k±Pp — 1, and 5vp is the rms 
amplitude of AW/KAW velocity fluctuations at A:xPp ^ 1. In 
writing Equation ( |47] i. we have taken the total fluctuation en- 
ergy per unit volume at ^xPp ^ 1 to be twice the kinetic en- 
ergy density p5vp/2, and we have assumed that dissipation 
at kj^ < Pp ' does not reduce the cascade power at kj^pp ^ 1 
much below the level that is present throughout the inertial 
range. There is some evidence that the magnetic fluctua- 
tions in the solar wind are co nsistent with turbulence theo- 
ries based on critical bala nce (iHorbury et al.l 120081: iPodestal 
l2009t iForman et aP 1201 ll) . However, there are conflicting 
claims in the literature over the validity of Equatio n (l46T l when 
^rms ^ ^rms' ^ point to which we return in Section l54l 

We assume that for length scales X between Pp and the per- 
pendicular AW correlation length (outer scale) L±, the rms 
amplitude of the AW velocity fluctuations at perpendicular 
scale A , is g<: X^^^ as in obs ervations at r = 1 AU (Podesta et al. 
I2007t IChen et aP 1201 ll) . direct numerical simulations of 
AW turbulence in the presence of a strong back ground 
magnetic field (jM ai'on & Goldreich 2001; Miiller & G rappinI 
120051; lMasonetal.1 i2006i: iPerez & Boldvrev 2008), and re- 
cent theories of strong AW turbulence feoldyrev. .2006: 
iPerez & Boldyrevll2009h . We flius take 



S,, — ^rms ( Pp 



1/4 



We set 



(48) 



(49) 



where Lxgj is the value of Lj_ at the coronal base, so that Lj_ 
increases in proportion to the cross-sectional radius of the flux 
tube in our model. 

In addition to linear damping, KAWs at k±Pp ^ 1 
undergo nonline ar damping through the "stochastic heat- 
ing" of protons (iMcCh esnev et al.1 Il987l; IChen et al.1 l200Tt 
iJohnson & Chen g 2001). In stochastic proton heating, 
AW/KAW fluctuations at k±Pp ^ 1 cause proton orbits in the 
plane perpendicular to Bq to become stochastic, and the pro- 
tons are subsequently energized by the time-varying electro- 
static potential. Using numerical simulations of test parti- 
cles int eracting with a spectrum of randomly phased AWs and 
KAWs. IChandran et alT(l2010h found that the effective damp- 



ing rate of KAWs at ^xPp ' 
proton heating is 



1 and piip < 1 from stochastic 



0.18epnp exp 



C2 



(50) 



where C2 is a dimensionless constant, 

_ Sv'p 
£p — , 
vxp 



(51) 



and vxn 
speed 



I p/fflp is th e proton perpendicular thermal 
Chandran etal] llOlO) found that C2 = 0.34 for ran- 



domly phased AWs and KAWs, but conjectured that C2 is 
smaller (and hence stochastic heating is more effective) in 
strong AW/KAW turbulence, because much of the dissipation 
in strong AW/KAW turbulence occurs within coherent struc- 
tures in which the fluctuation amplitudes are larger than their 
rms values (Dmitruk et al. 2004). In the vicinity of such struc- 
tures, proton orbits are more stochastic tha n on average, al- 
lowing for more efficient stochastic heating. FChandranl (120 1 01) 
developed a model of ion temperatures in coronal holes based 
on stochastic heating, and this model matched the observed 
temperature profile when C2 was set equal to 0.15. In our 
model, we leave C2 as a free parameter 
The total effective damping rate of KAWs at ^xpp = 1 is 



Ytot = ye+Yp + Ys. 



(52) 



We define T to be the fraction of the cascade power that is 
dissipated at fcxPp ^ 1- This fraction is roughly ytotfc when 
Jtotk ^ 1 and roughly 1 when ytotfc ^1- To interpolate 
smoothly between these limits, we set 



Ytotfc 



1 + Ytotfc 



(53) 



Landau damping and transit-time dampi ng of KAW s on pro- 
tons conti-ibute to g||p but not to Qxp (|Stixlll992h . On the 
other hand, stochastic heating leads primarily to perpendic- 
ular proton heating when p||p ^ 1 (Chandran et al. 2010). 
IJohnson & Chend (120011) have shown that stochastic heat- 
ing leads to significant perpendicular proton heating even 



at Pii 



1, and for simplicity we take stochastic heating to 



contribute only to Q±p, regardless of the value of P||p. We 
divide the cascade power that is dissipated at k±pp ^ I into 
three parts — electron heating, perpendicular proton heating, 
and parallel proton heating — in proportion to the correspond- 
ing damping rates. Ye, Ys, and Yp. As mentioned previously, we 
assume that the fraction of the cascade power that is not dis- 
sipated at A:xPp ^ 1 cascades to scales ^ Pp, dissipates via 
interactions with electrons, and contributes to Qe- This proce- 
dure leads to the relations 



(1+Yefc)e 
1+ Ytotfc ' 

Ysfcfi 



and 



^^P"T + Ytotf< 



Q\\v = 



1+ Ytotfc 



(54) 
(55) 

(56) 



3. NUMERICAL METHOD 

We integrate Equations ( fTOl i through ( fTTb forward in time 
until a steady state is reached using the implicit numerical 
method described by IHu et al.l (119971) . We use a logarithmic 
grid in r with grid points r, extending from 1 Rq to 1 .2 AU, 
where ; = 0, 1,2, . . . ,A^-|- 1 and = 1320. At each time step, 
we update the variables by integrating the equations only at 
ri,r2, . . . ,r7v. We then update the variables at ro and r^+i 



1 10 100 1 10 100 1 10 100 

r/Re r/R^ r/R^ 

Fig .2. — Left panel: Proton (or electron) number density in model (solid line). The x s show observed values in a polar coronal hole near solar minimum from 
Table 14.19 of Allen (1973). The filled circle is the mean proton density « = 2.7 ±0.86 cm"^^ at heliographic latitudes > 36° measured during Ulysses' first polar 
orbit, scaled to r = 1 AU (McComas et al. 2000). Middle panel: Important time scales in our steady-state numerical solution. Right panel: The ratios of thermal 
to magnetic pressure: p||p = 8ji«/:b7]|p/B^, pxp = ^T^nk^T^^^ / B- , and pe = Smk^T^/B'. 



using the following boundary conditions. At r = ro we set 
n = riQ, Te = Txp = 7||p = 7©, and £w = n0OTp(5v0)^, where 
Hq, Tq, and Sv© are constants (see Table [T]i. We determine 
U, q±p, and ^||p at r = ro by linearly extrapolating from the 
values at r = ri and r — r2- We determine C/(ro) in this 
way rather than by fixing the value of U (ro) so that the vari- 
ables can evolve towards a steady-state transonic solution that 
passes smoothly through the sonic point. We determine ^xp 
and ^iip at r = ro by linear extrapolation for the following rea- 
son. At radii slightly greater than ro, the values of ^xp and 
^IIp are determined to a good approximation by neglecting the 
terms on the left-hand sides of Equations (flSl l and ( fTSI l. be- 
cause the collision time Vp ' is much shorter than the expan- 
sion time r/U. We call the values of ^xp and ^||p determined 
in this way the "collisional values." If boundary conditions 
were imposed on ^xp and q^^^ at r — ro that were different 
from the collisional values, then in steady state a boundary 
layer would develop at r = ro of thickness ^ f/o/Vpo, where 
Uo = U{ro) and Vpo = Vp(ro), and q±^p and <7||p would ap- 
proach their collisional values at r ~ ro + t/o/Vpo- However, 
by linearly extrapolating the values of q±p and ^||p to r = ro, 
we prevent such unphysical boundary layers from appearing. 
At r = rA?xi we evaluate all variables by linearly extrapolating 
from their values at r = r^- i an d r — r^. For the solutions 
presented in Sections |4] and 15.71 we used the following ini- 
tial conditions: = Txp = 7||p = 7^,(3 - 2RQ/r){r/RQ)-^/'' , 
U = (655 km/s)[l +20{RQ/r)^]-\ n = nQUoaQ/{Ua), Ew = 
nOTp(5v0)^, and ^xp = q\\p — 0. 

4. A NUMERICAL EXAMPLE RESEMBLING THE FAST SOLAR 
WIND 

In this section, we focus on a steady-state solution to Equa- 
tions ( [Tol l through ([TtI i in which the model parameters are 
set equal to the values listed in Table [T] Our choice for 
L±Q is motivated by Faraday-rotation measurements along 
lines of sight passing through the corona and near-Sun solar 
wind0 (iHoUweg et al.ll20Tol) . After choosing the super-radial 

^ IHollweg et al.H2010l) found that the magnetic fluctuations SB in the solar- 
wind model of lCranmer et al] j2007l) led to close agreement with the fluctua- 
tions in the Faraday rotation of radio transmissions from Helios near superior 



expansion factors /max and Ri, v/e determine Bq from Equa- 
tion (|5]l and the condition 

B(l AU) =2.83nT, (57) 

which is the mean radial magnetic field strength mea- 
sured during Ulysses' first polar orbit, scaled to r = 
1 AU d'McComas et al.i2000l) . This yields = 11 .8 Gauss. 

TABLE 1 

Parameters in Numerical Example 



Quantity Value 

n@ 10* cm"3 

Tq 7 X 10^ K 

5v0 41.4 km/s 

/max 9 

Ri 1.29R0 

Lx0 10^ km 

Cd 0.75 

C2 0.17 

Hh 0.75 



The resulting numerical solution resembles the fast solar 
wind in several respects and illustrates the physical processes 
operating in our model. The density profile is shown in the left 
panel of Figure |2l Near the Sun, n decreases rapidly with in- 
creasing r. At larger r, as the wind approaches its asymptotic 
speed, n becomes roughly proportional to r^^. The Coulomb 
collision time scales v^t', Vpp'^;, and Vpg^ are shown in the mid- 
dle panel of Figure |2] Also plotted are the advection time 

conjunction. A key parameter in the model of Cranmer et al. (2007) was Lx, 
which was set equal to 75[(1470 Gauss)/B]'/2. For B = 11.8 Gauss (the 
value of B,., in our numerical solutions), this leads to = 837 km, which 
we have rounded to 10^ km in choosing the value of Lx, . ,. It is possible that 
Lx -i is significantly larger than this value, but in this case SB would have to 
be significantly smaller than in the model of Cranmer et al. (2007) in order to 
be consistent with the Faraday rotation measurements, preventing AW turbu- 
lence from providing the heating needed to power the solar wind. We exclude 
the possibility that L xa <jC 10^ km, b ecause then SB would have to be much 
larger than in work of lCranmer et al.l ^007) in order to be consistent with the 
Faraday rotation fluctuations, causing F„ to be much greater than the total 
energy flux of the solar wind. 
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scale fadv — r/U and the total proton scattering time scale Vp ' , 
which includes the effects of temperature-anisotropy instabil- 
ities. Close to the coronal base, the density is sufficiently 
large that Coulomb collisions play an important role, acting 
to maintain Tg ~ Txp ~ 7j|p. Farther from the Sun, however. 
Coulomb collisions cause only a negligible amount of proton- 
electron energy exchange and proton temperature isotropiza- 
tion. In the right panel of Figure |2] we plot the ratios of ther- 
mal to magnetic pressure, pe, Pxp, and p||p (defined in the 
figure caption), which vary from values <C 1 near the Sun to 
values larger than 1 at 1 AU. 

The profiles of the wind speed U and Alfven speed va are 
shown in the left panel of Figure[3] along with the rms ampli- 
tude of the fluctuating AW velocity 

5v™, = sj^. (58) 

In writing Equation dSST l, we have ignored the energy associ- 
ated with AWs propagating towards the Sun in the solar-wind 
frame, so that 5 

^rms — ^rms/^' our numerical solution, the 
Alfven critical point occurs at 

rA = 11.7/?o, (59) 

t/(rA) = 598 km/s, and U{1 AU) 800 km/s. The AW fluc- 
tuations lead to a turbulent heating rate per unit mass (2/p that 
is shown in the middle panel of Figure |3] 

The partitioning of the turbulent heating power between Qe, 
2xp, and 2||p is shown in the right panel of Figure[3] At r < 
\35Rq, electrons absorb most of the dissipated AW/KAW 
energy because 5vnns and hence gp are comparatively small 
(making stochastic heating weak) and because p||p ^ 1 (mak- 
ing proton Landau damping and proton transit-time damping 
weak). Between r = 1.35Rq and r = 1857^0, stochastic heat- 
ing is the most efficient dissipation mechanism in the model, 
and gxp > 0.5g. At r > 185/;0, Qe is the largest of the three 
heating rates. The parallel proton heating rate remains < 
at large r despite the fact that p||p > 1 and Yp > Ye. The reason 
for this is that at large r a significant fraction of the turbulent 
energy cascades to perpendicular scales <C pp at which the 
fluctuations dissipate on the electrons. 

In Figure |4] we plot T^, T±p, and 7||p, and in Figure |5] we 
plot the heating rates that help determine these temperature 
profiles. Between Rq and 2Rq, there are ripples in the plots 
of Qe and Q, which result from local flattenings in the Alfven 
speed profile (see Equations (|42] | and (|43]l)FI The ripples in 
Qc cause and to vary in such a way that the electron 
conductive heating rate partially offsets the variations in Qe- 
This can be seen in left panel of Figure |5] in which the elec- 
tron conductive heating rate — V • qe is positive at r > 2.25/?0 
but alternates sign at each sharp dip in the plot of | V • qe | • At 
Rq < r < 2.5/?0, AW turbulence is the dominant heat source 
for electrons, as the electron heat flux acts primarily to cool 
the electrons in this region. The left panel of Figure |5] indi- 
cates that a significant fraction of the turbulent heating power 
deposited into the electrons between Rq and I.SRq, is con- 
ducted away to either larger or smaller radii. At r > 3Rq, Qe 

' The oscillations in Q and Qe between r = S ;, and r ~ 2R result in part 
from our approximating z^,, and Q in Equations (42) and j43) based on the 
local value of \dv/^/dr\ . A more realistic treatment would account for the fact 
that AWs propagate some distance along the magnetic field before their 
energy cascades and dissipates, so that the local values of and Q depend 
upon the value of \dv/^/dr\ throughout some range of radii. 



and I V • qe| are of similar magnitude, although | V • qe|/2e 
grows to a value ~ 2 as r increases to 1 .2 A U. At r > IQQRp , 
the electrons approach the state described bv^ Hollwea (119761) . 
in which collisionless heat flux is the only source of electron 
heating and 

7;oc„2/[3(l+aH)l^ (60) 

However, the 7^ profile remains slightly flatter than the scaling 
in Equation (|60] | because of turbulent heating. 

At 2/?0 < r < IIRq, T±p is determined by a balance be- 
tween turbulent heating and solar-wind expansion, with the 
proton heat flux and collisions playing only a minor role, as 
shown in the middle panels of Figures |2] and |5] The quan- 
tity H±p plotted in Figure |5] is the perpendicular heating rate 
resulting from the proton heat flux, 

^xp = -^|:(«Vp)- (61) 

The Tj_p profile at 2/?0 < r < IIR,^ is also affected by the 
self-limiting nature of stochastic heating. As Tip increases, 
the fluctuations in the electrostatic potential energy at the 
proton-gyroradius scale become a smaller fraction of the av- 
erage perpendicular kinetic energy per proton, ^b7±p- As a 
consequence, these fluctuations have less effect on the pro- 
ton gyro-motion, and the proton orbits become less stochastic. 
This leads to a strong reduction in the stochastic heating rate 
when Tj_p exceeds a certain threshold that depends upo n 5vp 
and C2, as described in more detail bv I ChandranI (120 10 ). At 
r = IIRq, the plasma encounters the threshold of the oblique 
firehose instability, causing Vp to increase abruptly and lead- 
ing to sharp cusps in the T±p and 7||p profiles plotted in Fig- 
ured At r > IIRq the temperature anisotropy ratio 7xp/7||p 
evolves approximately along the oblique-firehose instability 
threshold, as shown in the right panel of Figure |4] 

The parallel proton heating rate associated with the proton 
heat flux, 

% = --rNip) + ^r' (62) 

adr a or 

is shown in the right panel of Figure |5] along with 2||p- Al- 
though |//||p| is small compared to Q, it is larger than 2||p at 
r < IORq. Despite the fact that |//||p| < Q, the proton heat 
flux causes 7]|p to increase with increasing r within this range 
of radii. This is possible because solar-wind expansion has 
only a small effect on 7||p at these radii. As can be seen from 

Equation ( fT4b . solar- wind expansion acts to make 7||p °^ ii^/B^ 
in the absence of competing effects. As the solar wind ap- 
proaches its asymptotic speed, n becomes approximately pro- 
portional to r^2. Likewise, when r exceeds a few Rq, B °^ r^^, 
at least in our model in which solar rotation is neglected (see 
Section lSTI i. When both n and B are °^ r^^, (double) adiabatic 
expansion neither increases nor decreases 7||p. 

The electron Coulomb mean free path A-mfp and electron 
heat fluxes qe are plotted in the left and middle panels of Fig- 
ure |6l The electron heat flux transitions from the collisional 
regime to the collisionless regime in our model at r = rn = 
5Rq, approximately the point at which A-ipfp = r/2 as in the 
coUisionless-heat-flux model of Ho llwed fl974ll976l) . In the 
collisionless regime, qe is smaller than the Spitzer-Harm heat 
flux qe s and somewhat smaller than the free-streaming heat 
flux 

^sat.e = l.SnkBTeVte, (63) 



9 




Fig. 3 . — Left panel: Th e open box represents the range of upper hmits on 5i'rn,s obtained bv lEsser et al.' fl999^. The filled circles are Helios measurements 
of8\'nns <Bavassano eT aP 2000). The filled square is the median proton flow speed U = 761 km/s at heliographic latitudes > 36° during Ulysses' first polai' 
orbit, scaled to r = 1 AU iMcComas et al. 2000). The size of this square represents the range 702 km/s < U < 803 km/s corresponding to the 5* through 95* 
percentiles of the measured distribution. Middle panel: Total turbulent heating rate per unit mass. Right panel: The fractions of the turbulent heating power that 
go to electron heating, perpendicular proton heating, and parallel proton heating. 
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Fig. 4. — Left panel: Triangles are electron temperatures inferred from spectroscopic observations of a polar coronal hole ('Landi'20080. The square is the mean 
electron temperature in ISEE 3 and Ulysses measurements of fast-wind streams with 600 km /s <U < 70 km/s iNewburv et al. 1998). Middle panel: Triangles 
are proton kinetic temperatures inferred from UVCS measurements of a polar coronal hole jEsser et aT. 1999). Open squares (xs) are perpendicular (parallel) 
proton temperatures measured by Helios in fast solar-wind streams with 700 km/s <U < 800 km/s (Marsch et al. 1982b). Right panel: The value of Tj_j,/Ty in 
our numerical solution progresses from left to right along the solid-line curve as r increas es from l/?,.;, to 258^0 . The dotted lines give the instability thresholds 
for the mirror instability (upper curve) and oblique firehose instability (lower curve) from lHellinger et all <2006[) . at which the maximum instability growth rates 
equal lO^'flp, where I2p is the proton cyclotron frequency. The long-dashed line is the fit 71p/7]|p = l-l^Pyp"^'^ obtained bv iMarsch et alj t2004l) to Helios data 



of high-speed wind streams between r = 0.29 AU and r = 0.98 AU. 

but comparable to Helios measurements of the ele ctron heat 
flux in the fast solai" wind jMarsch & Richteiill984l) . 

The proton heat fluxes are plotted in the right panel of Fig- 
ure |6] As this figure shows, ^xp ™d q^^^ are significantly 
smaller than the free-streaming heat flux 



^sat.p = l.SnksTpVtp, 



(64) 



where Vtp = V^fiTp/'Mp- This can be understood on a qual- 
itative level from the following argument. If a proton tem- 
perature gradient were set up in a collisionless plasma with 
no background flow and no initial heat flux, then the proton 
heat flux would grow in time, approaching a level compara- 
ble to the free-streaming value after a time ^ fcross — h/ ^tp, 
where Ij = Tp/j VT^pj. In our model, Ij ^ r, and this "cross- 
ing time scale" is a factor of ^ M larger than the expansion 
time scale of the solar wind, fadv — r/U, where M = [//vtp is 
the Mach number. In our numerical solution, M equals 3.56 
at r = IQRq and grows monotonically with increasing r to a 



value of 15.1 at r = 215Rq. Thus, throughout most of our 
solution, fadv — r/U fcross- As a result, the protons in our 
model do not have time to set up a heat flux comparable to 
the free-streaming heat flux within the time it takes for the 
plasma to double its distance from the Sun, which reduces 
^Xp ™d ^iip relative to their values in a stationary plasma with 
comparable density and temperature profiles. 

As discussed in Section |2] the total energy flowing through 
the flux tube in our model per unit time, flFtot, is independent 
of r in steady state. (In our numerical solution, the ratio of the 
maximum to minimum values of aFtot is 1.003.) This makes 
it straightforward to identify the principal source of energy in 
our model and to understand how energy is converted from 
one form to another as plasma flows away from the Sun. The 
total energy flux defined in Equation (l27T i is the sum of the 
bulk-flow kinetic energy flux 



(65) 
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Fig. 5. — Heating rates associated with the heat fluxes and the dissipation of AW/KAW turbulence. 




Fig. 6. — Left panel: A.„,fp is the electron Coulomb mean free path. Middle panel: The electron heat flux q^, free-streaming heat flux ^^sat.o and Spitzer-Harm 
heat flux q^s,- The circles are Helios measurements of the electron heat flux in high-speed wind with U > 600 km/s i Marsch & Richter 1984) . Right panel: The 
proton heat fluxes (yxp and (/[p and free-streaming heat flux qsM,p. 



the gravitational potential energy flux 

UGMqP 



the enthalpy flux 



UnkB 



5n 



37||, 



the total heat flux 



qtot^qe + q±p + q\\p 

and the AW enthalpy flux 



(66) 



(67) 



(68) 



(69) 



We plot these fluxes in Figure [T] normalized to the total en- 
ergy flux Flat- As this figure shows, the wind in this solution is 
driven fundamentally by the AW enthalpy flux. As the plasma 
flows away from the Sun, part of the AW enthalpy flux is con- 
verted into gravitational potential energy flux as the flow lifts 
material out of the Sun's gravitational potential well. Most of 



the remaining AW enthalpy flux is gradually converted into 
bulk-flow kinetic energy flux, which dominates the total en- 
ergy flux at r = 1 AU. 

5. DISCUSSION 

In this section, we discuss our results and present a second 
steady-state solution that incorporates pitch-angle scattering 
by the cyclotron instability. 

5.1. Solar Rotation 

If solar rotation were taken into account, B would follow 
the Parker spiral, and at large distances from the axis of rota- 
tion B would become approximately azimuthal rather than ra- 
dial, withB o<: r^' instead of B 0= r^^. Assuming n °^ r^^, dou- 
ble adiabatic expansion in this azimuthal-field regime would 
lead to the scalings Txp °^ r^^, ^|p °^ ^'^'^ 7xp/7||p °= r. 
In contrast, in the radial magnetic field of our model, dou- 
ble adiabatic expansion leads to Txp °^ r^^, 7||p °^ r", and 
rxp/7||p 0= r^^ assuming n °^ r^^. The inclusion of rotation 
would thus increase the temperature anisotropy ratio Txp/T^ip 
at large r. In the solar wind, the transition between radial 
and azimuthal magnetic field occurs gradually throughout a 
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Fig. 7. — The fractions of the total energy flux F,ot that come from the 
bulk-flow kinetic energy flux Fu, the AW enthalpy flux F„. the gravitational 
potential energy flux Fg, the total heat flux qtat, and the enthalpy flux F^. 



range of radii centered at r± ^ U /Q., where Q. is the angular 
frequency of the Sun's rotation, and rj_ is distance from the 
Sun's spin axis. For fast wind with U = 800 km/s, and for 
n = 2.64 X 10'^ s~^ (the v alue of H at a solar latitude of 45° 
(iSnoderass & Ulrich|[l990h '). U/Q.= IA AU. 

In addition to modifying the temperature anisotropy ratio 
at large r, the inclusion of rotation would increase the to- 
tal magnetic field strength at large r, thereby reducing P||p. 
For example, during the first polar orbit of Ulysses, the ra- 
tio of the mean total field strength (scaled to r = 1 AU) 
to the mean radial field strength (scaled to r = 1 AU) was 
~ 1.7 (McComas et al. 2000). At fixed « and 7||p, increas- 
ing B by a factor of ~ 1 .7 reduces p||p by a factor of ~ 3. Our 
overestimate of p||p at r 1 AU causes us to overestimate 2||p 
and also pushes the thresholds of the firehose and mirror in- 
stabilities towards smaller temperature anisotropies. 

5.2. The Need for Parallel Proton Cooling 

In the numerical solution presented in Section|4] 7||p > 7xp 
at r > 35Rq. In contrast, in Helios measurements of fast- 
wind streams with U > 700 km/s, T±p typically exceeds 7||p 
between r 60/?© and r = 130/?© (Marsch e t allll982b) . Part 
of this discrepancy may be due to our neglect of solar rota- 
tion (Section 15. Il l or one or m ore perpendicular proton heat- 
ing mechanisms (Section 15.6b . This discrepancy may also 
arise, at least in part, from our neglect of kinetic mecha- 
nisms that act to reduce 7j|p. For example, part of the par- 
allel proton thermal ener gy in the fast solar wind is in the 
form of a proton beam ( Marsch et al.lll982bl |2004 . When 
the relative speed of the beam component with respect to 
the core of the proton distribution exceeds ~ va, the pro- 
ton beam excites plasma instabilities that slow the beam 
down (Daughton & Gary 1998). Since va decreases with in- 
creasing r, these instabilities lead to the steady deceleration of 
the proton beam compon ent, which reduces the p arallel pro- 
ton thermal energy (Hell inger & Travmcekll20l"11) . The pos- 
sible need for parallel proton coohng in the solar wind was 
previously suggested bv lHu et aTl (119971) . 



5.3. Uncertainties in the Total Turbulent Heating Rate 

Two of the assumptions in our estimate of Q become in- 
creasingly inaccurate as r increases. First, we have assumed 
that ^ i.e., that most of the AWs propagate away 
from the Sun in the sola r-wind frame. While zty,^ is likely 
Z^c in c oronal hol es (ICranmer & van Ballegooiienll2()05l: 
IVerdini & V elli 2007; ICranmedl2010r Helios measurements 
in the fast solar wind show that Zrms/^rms decreases fro m :^ 4 
to_^2 as r increases from %ORp, to 1 AU (Bava ssa no et al.l 
120001) . Second, the model of lChandran & Hollwed (120091) that 
we employ to estimate assumes that the energy cas- 
cade time at the outer scale, t^^^^. ~ Lx/^rms' much shorter 
than the linear wave period P. For the numerical solution pre- 
sented in Section m t^^^,. grows steadily as r increases, reach- 
ing a value ~ 3 x lO^* s at r = \OORq and a value ~ 10"* s 
at r = 2OO/?0. Thus, for P ^ I hour, the assumption that 
&c ^ P breaks down at large r. We also note that veloc- 
ity shear may be an important additional source of AW tur- 
bulence in the^ solar wind, one that is not included in our 
model (Ro berts et al.lll987 ). As a source of both z+ and 
fluctuations, AW excitation by veloc ity shear acts to dec rease 
the ratio Zrl^s/^rms in solar wind ([Roberts et al.lll992l) . 

5.4. Uncertainties in the Division of the Turbulent Heating 
Power Between Qe, Q±p, and Q\\.g. 

Our prescription for partitioning the turbulent heating 
power between protons and electrons, and between paral- 
lel and perpendicular proton heating, depends strongly upon 
two of the assumptions we have made in modeling AW 
turbulence in the solar wind: the scaling of 5vp in Equa- 
tion ( |48] | and the "critical balance" condition in Equation (|46] |. 
Equation (|48] | corresponds to an assumption that the inertial- 
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range velocity power spectrum is, k^^ , wh ich is consis- 
tent with spacecraft measurements at r = 1 AU (Podest a et al.l 
2007; Podesta & Bhattacharjee 2010; Chen et al. 201 1|) and 
direct numerical simulations of AW turbulence in which the 
fluctuating magnetic field is < 0.2 times the backg round 
magnetic field jM aron & Goldreich 2001; Miiller & G rappinI 
120051; iMason et al. 2006; Perez & Boldvrevi i2008i) . includ- 
ing simulati ons t hat account for cross helicity, in which 
~rnis > ^rms dPerez & BoldvrevI 120091) . The power spec- 
trum of the magnetic field is typically steeper (~ k^^l^) 
than th e velocity power spectrum in spacecraft measure- 
ments (Matthaeus & Goldstein '1982'; 'Goldstein et al.' '19951 
Bruno & Carbone 2005; Podesta etal. 2007; Chenet^ 
2011) and direct num erical simulations dMiiller & Grappml 
,2005>:iBoldyrev et al.l201 1 ). The stochastic heating rate in our 
model, however, depends upon the velocity power spectrum, 
not the magnetic power spectrum, because the velocity spec- 
trum is a measure of the electric-field fluctuati ons, which con- 
trol the stochastic heating rate when p||p < 1 (IChandran et al.l 
i2010) . On the other hand, AW turbulence near the Sun is 
strongly affected by non-WKB wave reflection and may dif- 
fer from AW turbulence at r ~ 1 AU and from the homo- 
geneo us turbulence in the numerical simulations mentioned 
above ( IVerdini et al.l2009h . If the velocity spectrum is steeper 
(flatter) than we have assumed, then the amplitude of the ve- 
locity fluctuations at ^xPp ^ 1 is smaller (larger) than in Equa- 
tion ( |48] |, and stochastic heating is weaker (stronger) than in 
our model. 

The critical-balance condition that we have adopted in 
Equation (|46] | is the same as the condition in Boldyrev's 
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(2006) theory of scale-dependent dynamic alignment and is 
consistent with Perez & Boldyrev's (2009) extension of this 
theory to the cross-helical case, which holds that the cas- 
cade times of z+ and fluctuations are the same even when 
Zrins-^Zms- Equation (|46] | is also consistent with the work 
of iPodesta & Bhattacharjee (2010) if the ratio p/q in their 
analysis is set equal to unity. However, Equation ( l46b is not 
consistent with three competing (and mutually incon sistent) 
models of AW turbulence wi th cross helicity ( Lithwic k et alj 
I2007t [Beresnyak & Lazai-ian"2008^, IChandranll2008ai) . More- 
over, none of the six studies just mentioned accounted for the 
non-WKB reflection of AWs, which could affect wavenumber 
anisotropy in the solar wind. The range of parallel wavenum- 
bers that are present at ^xPp = 1 in AW turbulence in the solar 
wind is thus uncertain. If the typical values of |^|| | are larger 
than we have assumed, then linear wave damping is stronger 
than we have assumed, implying that 2±p/2 is smaller than 
in our model. In addition, near r = 1 AU, increasing the lin- 
ear damping rates would increase Q\\^/Qe, because less power 
would cascade to scales ^ Pp, and because protons absorb 
most of the power that is dissipated at ^xPp ~' 1 when pp > 1. 

Other sources of uncertainty in our partitioning of the turbu- 
lent heating power include the stochastic-heating constant C2 
in Equation ( fSOl l. which is not known for the case of pro- 
tons interacting with strong AW/KAW turbulence. If we 
have underestimated C2, then we have overestimated the ef- 
ficiency of stochastic heating, causing our model to overesti- 
mate Txp. In addition, the use of linear Vlasov theory to esti- 
mate damping rates may be inaccurate when applied to large- 
amplitude AW/KAW turbulence in the solar wind (Borovsky 
& Gary 20 11; but see also Lehe et al. 2009). 

5.5. The Coronal Electron Temperature 
Electron temperatures in the low corona inferred from 
line ratios are found to be in the range ~ 8 x 10^ - 
10*^ K (Habbal et alJ [T99l IPosch ek et al. l 200Tt IWilhelml 
120061: iLandii l2008h . similar to the electron temperature in 
our model. On the other hand, Ulysses measurements of 
ion charge states in the fast solar wind emanating from the 
south polar coronal hole suggest that T^ reaches a maximum 
of ~ 1.5 X 10*^ K at r - 1.3 - \.5Rq (IKo et al.i ri997). The 
numerical solution we have presented in Section |4] does not 
reach such high electron temperatures, which may indicate 
that we have under-estimated the electron heating rate. On 
the other hand, the electron temperatures inferred from ion 
charge-state ratios at 1 AU may be inflated to some degree by 
the presence of superthermal electrons in the corona, which 
could significantly en hance the rates of ionization into high- 
energy charge states (lOwocki & Scudded [19831 lBurgilll987l: 
iKo et al.lll996h . potentially making the T^ profile in the left 
panel of Figure |4] consistent with the ion-charge-state mea- 
surements. 

5.6. Other Possible Heating Mechanisms 

Although low-frequency AW turbulence is the only non- 
conductive heating mechanism in our model, other mecha- 
nisms may be important in the solar wind. For example, 
compressive waves are believed to play an important role in 
chromospheric heating, and may deposit a significant amount 
of energy in the l ow corona as well (iCranmer et al.l l2007t 
IVerdini et al.ll2010l) . Type-II spicules may also be an impor- 
tant s ource of heating in the low corona (i De Pontieu et al.l 
120111) . Farther from the Sun, the solar wind may be heated 



by high-frequency waves that are generated by either a turbu- 
lent cascade ([L eamon et al. 1998; Hollweg & Isenberg 200!^ 
[Hamilton et al.ll2008i; Jsenberg & Vasquezii201 L) or by insta- 
bilities driven by the differential flow between the core of the 
proton velocity distribution and either alpha particles or pro- 
ton beams ( Gary et al. 2000; Hellinger & Tray nfcek 201 1). 

5.7. Cyclotron Instability versus Mirror Instability 

iHellinger et al.l (l2006h showed that near r = 1 AU , the 
proton temperature anisotropy ratio 7xp/7||p is limited from 
above by the threshold of the mirror instabil ity. However, us- 
ing Helios measurements, Bou rouaine et al.l (12010) found that 
rxp/7||p appears to be limited from above by the cyclotron 
instability in fast solar wind streams between 0.3 AU and 
0.4 AU. To investigate how this latter possibility would af- 
fect our numerical solutions, we have repeated the numerical 
calculation presented in Section |4] with the parameter values 
listed in Table [T| using a model for enhanced proton pitch- 
angle scattering based on the cyclotron instability instead of 
the mirror insta bility. That is, we hav e set R^i Rc in Equa- 
tion (O, with (.HelUnger et al.ll2006l) 

= 1 +0.43(P||p + 0.0004)"''-^2. (70) 

The density and temperature profiles in this second steady- 
state solution are shown in Figure |8] The density profile is 
very similar to the profile in the left panel of Figure |2] The 
profiles of Tf. and Txp are very similar to the profiles shown in 
the left and middle panels of Figure |4] The most notable dif- 
ference between the temperature profiles in the two numerical 
solutions is in 7||p, which is larger at 4Rq < r < IQR^ when 
the cyclotron instability threshold is used. As can be seen 
in the right panel of Figure [8] when 7xp/7||p is limited from 
above by the anisotropy ratio in Equation (iTOl i. the plasma en- 
counters and then evolves approximately along the cyclotron- 
instability threshold at small r and small p||p. At larger r and 
larger P||p, the plasma evolves approximately along the thresh- 
old of the oblique firehose instability, as in the numerical so- 
lution presented in Section|4] 

5.8. Sharp Boundaries in Parameter Space 

We have found that there are regions in parameter space 
in which tiny changes in certain parameters lead to large 
changes in the final steady-state solutions. This phenomenon 
may be related to the abrupt transition from fast wind to 
slow wind in theoretical models in which t he magnetic ge- 
ometry varies with heliographic latitude (ICranmerl 120051: 
iCranmer et a l. 2007). 

5.9. Comparison to Previous Studies 

A number of authors have developed s olar-wind mod- 
els incorporati ng t emperature anisotropy (Leer & Axford 
[1971 IWha ng..l97l lDemars & Schunkl[T99Tl Hu et al. 199% 
Olsen & L eerll999i;lEndeye & Leerf2001l;lLie- Svendsen et at] 
1200 U iJanse et al.1 l2006h . and some of these models 
also included energy and/or momentum deposition by 
AWs ( Uu et al.l[T997t l01sen & Leeill999tlLie^vendsen et al.1 
2001). Our work, however, has several features that are not 
present in these previous studies. First, we evaluate the total 
turbulent heating rate Q using an analytical theory of low- 
freq uency AW turbulence driven by non-WKB wave reflec- 
tion dPmitruk et al.ll2002l: IChandran & Hollweg||2009h . Sec- 
ond, we divide the total turbulent heating rate Q into three 



13 




1 10 100 1 10 100 -4 -2 



r/Re r/R^ logio(^?||) 

Fig. 8. — Steady-state numerical solution when r^p/7]|p is limited from above by the threshold of the cyclotron instability (Equation ilOi ) rather than the mirror 
instability. The parameter values for this solution are given in Table[T] Left panel: Same as left panel of Figure|2] Middle panel: Symbols have same meaning as 
in the middle panel of Figure|4] In addition, the circles are electron temperatures infen'ed from spectroscopic observations of a polai' coronal hole I Landi 200^ 
and th e filled square is the mean electron temperature in ISEE 3 and Ulysses measurements of fast-wind streams with 600 km/s <U < 700 km/s (Newbury et aU 
[T99l) . Right panel: Same as right panel of Figure|4] except that the upper dotted-line curve is the instability threshold for the cyclotron instability in Equation ITOt T 



parts (Qe, 2xp ™d 2||p) using an analytical model of the 
collisionless dissipation of low-frequency AW turbulence. 
Third, we do not include any heating other than that which 
is provided by the proton and electron heat fluxes and low- 
frequency AW turbulence. Fourth, we account for the mirror 
(or cyclotron) and oblique firehose instabilities by enhancing 
the proton pitch-angle scattering rate when the proton temper- 
ature anisotropy exceeds the threshold of either instability. 

Two other groups of authors have developed solar-wind 
models based on energy and momentum deposition by low- 
frequency AW turbulence driven by non-WKB wave reflec- 
tion (Cranmer et al. 2007; Verdini et al. 2010). These au- 
thors allowed for a broad spectrum of AW frequencies at the 
coronal base and accounted for the fact that higher-frequency 
waves undergo less reflection. In contrast, our estimates for 
and Q in Equations (|42] | and ( |43] | were obtained in the 
low-frequency limit (Chandran & Hollweg 2009). On the 
other hand, these models treated the solar wind as a sin- 
gle fluid, without distinguishing between proton and electron 
temperatures. Also, neither model incorporated temperature 
anisotropy. 

6. CONCLUSION 

We have developed a two-fluid model of the solar wind that 
accounts for proton temperature anisotropy, pitch-angle scat- 
tering from mirror (or cyclotron) and firehose instabilities, 
and (kinetic) Alfven Wave (AW/KAW) turbulence. We ne- 
glect rotation and consider solar wind flowing along a narrow 
magnetic flux tube centered on a radial magnetic field line. 
The turbulent heating in our model is partitioned between 
electrons and protons, and between perpendicular and parallel 
proton heating, in accord with recent results on stochastic ion 
heating (Chand ran et al.l l2010) and linear damping rates cal- 
culated from the full hot-plasma dispersion relation. The elec- 
tron heat flux in our model transitions from the Spitzer-Harm 
value in the collisional region near the Sun to the Hollweg 
( 119761) value in the nearly collisionless conditions at larger r. 
To evaluate the proton heat flux in the presence of temperature 
anisotropy in both collisional and collisionless conditions, 
we use a fourth-velocity-moment fluid closure based on the 



guidi ng-center Vlasov equation (iKulsrudll 19831; ISnvder et al.l 
119971) . Our model conserves energy, and in steady state the 
total energy flowing through the cross section of the flux tube 
becomes independent of r. No energy is added to the solution 
through ad hoc heating terms. 

In Section|4]we present a steady-state numerical solution to 
our model equations for Rq < r < 1.2 AU, which we analyze 
in detail in order to gain insight into the different physical pro- 
cesses operating within the model. As shown in Figure [T] it 
is the AW enthalpy flux F^, that drives the solar wind in this 
solution. As plasma flows away from the Sun, the AW en- 
thalpy flux is gradually converted into gravitational potential 
energy flux and bulk-flow kinetic energy flux. By the time the 
flow reaches 1 AU, the total energy flux is dominated by the 
bulk-flow kinetic energy. As shown in Figure |5] electrons are 
heated primarily by the dissipation of AW/KAW turbulence 
at r < 2.5Rq. At r > 37?0, AW turbulence and conduction 
provide comparable amounts of electron heating, to within a 
factor of - 2. Between 1 .35Rq and 185/?©, AW/KAW turbu- 
lence dissipates primarily via stochastic proton heating, lead- 
ing to substantial perpendicular proton heating, as shown in 
the right panel of Figure |3] Parallel proton heating via Lan- 
dau and transit-time damping accounts for < 10% of the total 
turbulent heating at all radii in this solution, and g||p < Q.QIQ 
at r < lOO/?0. The proton heat fluxes ^xp and ^||p cause 7||p to 
increase gradually as r increases from 2Rq to ~ 207^0, despite 
the fact that solar-wind expansion reduces q±p and q\\p rela- 
tive to their free-streaming values. At r HR^s, the plasma 
encounters the threshold of the oblique firehose instability, 
causing the proton pitch-angle scattering rate to increase. At 
r > 71/?0, the proton temperature anisotropy ratio evolves ap- 
proximately along the firehose-instability threshold as p||p in- 
creases. 

This numerical solution is broadly consistent with a number 
of observations, as illustrated in Figures |2] |3] IH and |6] sup- 
porting the idea that AW turbulence may be one of the primary 
mechanisms responsib le for heating coronal holes and accel- 
erating the solar wind (iParkerll 1 9651: IColemanll 1968 ). Perhaps 
the most notable achievement of our model is that it comes 
close to explaining observations of perpendicular proton tem- 
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peratures, even though the heating in the model is provided by 
low-frequency AW turbulence rather than resonant cyclotron 
interactions. The main uncertainties in our results are asso- 
ciated with the stochastic heating rate in strong AW/KAW 
turbulence, the wavenumber anisotropy {k^'^/k^) in reflection- 
driven AW/KAW turbulence, the total turbulent heating rate 
at large r, and the effects of solar rotation on the temperature 
anisotropy ratio rxp/7||p. 

One of our principal objectives in this work has been to to 
connect theoretical studies of microphysical processes with 
observations of macrophysical quantities in the solar wind. 
By comparing our model to observations, we have been able 
to obtain a new test of the viability of AW turbulence as a 
mechanism for heating the solar wind and coronal holes. At 
this point, the results of this test are encouraging, but not 



fully conclusive, because of the uncertainties described above. 
However, as our understanding of kinetic plasma physics and 
turbulence in the solar wind progresses, it will be possible 
to use models such as the one we have developed to obtain 
increasingly rigorous tests of competing theories and, ulti- 
mately, to gain greater insight into coronal heating and the 
origin of the solar wind. 
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APPENDIX 



FOURTH-MOMENT FLUID CLOSURE OF THE GUIDING-CENTER VLASOV EQUATION 



The d erivation of th e equations in our model begins with Kulsrud's formulation of collisionless MHD for a proton-electron 
plasma ( lKulsrud ''1983'). Kulsrud's approach was to expand all quantities in the Vlasov and Maxwell equations in powers of 1/e, 
where e is the proton charge, and to consider the limit e ^ oa. This limit corresponds to the case in which the Debye length A-d 
and proton gyroradius Pp are much smaller than the length scales over which the macroscopic quantities vary appreciably. The 
fundamental variables in Kulsrud's theory are the mass density p, the fluid velocity U (which is the same for electrons and protons 
to lowest order in 1 /e), the magnetic field B, the proton and electron distribution functions /p and /e, and the parallel component 

of the electric field E, given by En =b - E, where b = B/B. To lowest order in 1/e, these variables satisfy the equations 



| + V.(pt/)^0, 



rdU \ (VxB)xB 

^ = Vx(C/xB), 
P = Y,pUI-bb)+J^PllM, 



(Al) 
(A2) 

(A3) 
(A4) 



f,\v-vub-VE\^d^v, 



f,{v\\-U\\fd\, 



^e.5«, = 0, 



and (see Snyder et al. 1997) 
d 



^^{m + V ■[f,B{v\\b + VE)] + ^ 



fsB{ -b -^-ljb VB 



Dt 



0, 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



where s is a subscript indicating particle species (p for proton and e for electron), is the distribution function of particle species s. 



nis and gj are the mass and charge of species s, v is particle velocity, vii = 6 • v, ve — c{E x B)/B , /j 



v\\b-VEY-/2B, 



= U ■ b, and D/Dt ~ d/dt + (v|| b + ve) ■ "V . In Equation iA% . fs is regarded as a function of position x, time t, magnetic 
moment /j, and parallel velocity vii . Rather than retain the subscripts on the number densities, we define 



(p, 



(AlO) 

which is also equal to because of Equation ( IA7b . We have neglected the electron contribution to the mass density, setting 
p = «mp. The parallel and perpendicular temperatures are related to the parallel and perpendicular pressures defined above in the 



usual way; pxs — «^b7±s and pn 
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J 
dGi 



Snyder et alj 



Oross & Krook 



1997h extended Kulsrud's collisionless M HD t o account for collisions by adding a BGK collision operator 
TT95I to the right-hand side of Equation iA9i . For the case we consider here, in which the electrons and 



ions have the same average velocity U and number density «, this collision operator takes the form 



C(/v)=Ev,,(Fm.-/v) 
k 



where 



\2%kBTs 



3/2 



exp 



ms(v|| - Ms/jB 



2kBT, 



kBTs 



is a shifted Maxwellian with temperature 



271, + 7||, 



(All) 



(A12) 



(A13) 



and Vsii is the collision frequency for momentum exchange bet ween species s and s pecies k. (Here, we neglect energy exchange 
between protons and electrons, but we include it in section |2]) ISnyder et alJ (Il997h then obtained a hierarchy of fluid equations 
by multiplying Equation iA9i by various powers of vy and fi and then integrating over V[| and /j. For the protons, the equations 
for /7xp and can be written jSnvder et alJI 19971 iSharma et al.ll2006l) 



Pfi^ (^) = - V • (?xpb) - ^xp V • b + ^ (Piip - /.xp) 



and 



where 



2B^dt 



3^ =-V-(<7||pb)+<?xpV-6+^(pxp-P||p), 



nip is the proton mass, and 



(A14) 
(A15) 

(A16) 
(A17) 

(A18) 



In Section|2]we neglect the Vpe term in Equation (IA18b because it is smaller than the proton-proton Coulomb collision frequency. 
Our ^iip is by definition a factor of 2 smaller than Snyder et al.'s (1997). In this appendix, the Lagrangian time derivative is given 
by 

| = | + [/.V, (A19) 



which generalizes Equation ^ to allow for arbitrary flow velocities. The equations for q±p and ^||p are 



and 



where 



dt \ p 



dt I p4 



PXp(P||p-/'Xp) 



-'"XX -nix 



-Vp^ll 



2 ^ II" ^ 2p 



3 

P+2 



P\\v{P\\p- Pi-v) 



Vb 



:p| /p^B(v||-f/||)2af^, 



and 



(A20) 
(A21) 

(A22) 
(A23) 

(A24) 



In Equation ( IA20b . we have corrected a minor error in Equation (19) of lSnyder et alJ (Il997h : the fourth term on the left-hand side 
of their Equation (19) should be qj_s^ ■ U instead of ^xsV • {Ui^\b). 
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To close these fluid equations, we set /p — Fbm in Equations ( IA22b through ( IA24b . where 



3/2 



Fbm — 



exp 



Mp/jB mp(v|| - f/||)2 



(A25) 



is a bi-Maxwellian distribution with the same values of «, U\\, Txp and 7]|p as the protons. This enables us to rewrite Equa- 
tions (IA20b and (IA21b as 



dt \ J 



-Vp^Xp ^ 



nkiT\\, 



-bVT 



-Lp- 



«^B^±p(7±p-7||p) 



Vb 



and 




3«feg7||j 
2m„ 



■bVTn, 



(A26) 



(A27) 



Equations (IA26b an d ( IA27b we re derived in d ifferent ways and with differing treatments of collisions bv lEndeve & Leej (1200 Ih . 
iLie-Svendsen et alj (12001 ), and lRamosl(l2003h . 

For the electrons, we set T±g = 7^|e = Te as described in section [21 and close the electron fluid equations by specifying the 
electron heat flux in terms of lower moments of the electron distribution (Section r2.1b . The electrons then satisfy a standard 
energy equation. 



2" 



5/3; 



(A28) 



In Section|2] we adapt the general equations given in this appendix to our ID solar-wind model and add extra terms to incorporate 
the effects of AW turbulence, collisional energy exchange between protons and electrons, and temperature isotropization by 
firehose and mirror instabilities. 
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